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ABSTRACT 

A model based on disk-stability criteria to determine the number of spiral arms of a general disk 
galaxy with an exponential disk, a bulge and a dark halo described by a Hernquist model is presented. 
The multifold rotational symmetry of the spiral structure can be evaluated analytically once the 
structural properties of a galaxy, such as the circular speed curve, and the disk surface brightness, 
are known. By changing the disk mass, these models are aimed at varying the critical length scale 
parameter of the disk and lead to a different spiral morphology in agreement with prior models. 
Previous studies based on the swing amplification and disk stability have been applied to constrain 
the mass-to-light ratio in disk galaxies. This formalism provides an analytic expression to estimate 
the number of arms expected by swing amplification making its application straight-forward to large 
surveys. It can be applied to predict the number of arms in the Milky Way as a function of radius and 
to constrain the mass-to-light ratio in disk galaxies for which photometric and kinematic measurements 
are available, like in the DiskMass survey. Hence, the halo contribution to the total mass in the inner 
parts of disk galaxies can be inferred in light of the ongoing and forthcoming surveys. 

Subject headings: Galaxy: disk - Galaxy: evolution - galaxies: kinematics and dynamics - stars: 
kinematics and dynamics 


1. INTRODUGTION 

A long-standing problem in understanding the dynam¬ 
ics of disk galaxies is the uncertainty of the fractional 
contribution of the dark halo to the total mass of the 
galaxy within the optical radius of the stellar disk. While 
galaxies are believed to be dark matter dominated in 
their outer parts, at radii much larger than the optical 
radius, the assessment of the relative fraction of baryons 
and dark ma tter in the inner part o f disk galaxies is still 
debated (e.g. |Gourteau et al.||2014|). 

The circular-speed curve constrains the total mass dis¬ 
tribution at all radii, and the disk contribution is usually 
estimated by combining the light profile with an estimate 
of the mass-to-light ratio of the disk. The difference be¬ 
tween the total rotation curve and the circular velocity 
of the d isk is assumed to be due to the presence o f dark 
matter (Bosma 19781 Carignan & Freeman 1985). Be¬ 
cause the mass-to-lignt ratio of the disk is uncertain, the 
relative contributions of the disk and halo are poorly de¬ 
termined. 

Attempts to constrain the mass-to-li ght ratio com e 
from studies of stellar populations (e.g. Conroy 2013), 
and comparisons bet ween dynamical and ste llar popnla- 
tion mass estimates (|Bell & de Jong||2001|). However, 
these models depend on assumptions made about the 
IMF and stellar ages, leading to rather large uncertain¬ 
ties in the disk contribution. A maximum disk hypothesis 
for the mass decomposition has b een suggested for disk 
galaxies (van Albada et al. 1985). It assumes that the 
total rotation curve of disk galaxies out to the optical 
radius can be fitted by models that assume a mass-to- 
light ratio of the disk independent of radius and with 
no need for dark matter in the inner parts (maximum 
disk). A typical maximum disk contributes more than 
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80% of the total circular speed at 2.2 R^, where R^^ is 
the scale length of the exponential disk. In terms of disk 
mass fraction the maximum disk hypothesis requires a 
disk contribution of /d(2.2 Rd)={ydisk/^tot)^ >0.7. 

Additional arguments have been made to constrain 
the halo and disk contribution to the total gravitational 
field via different techniques. In particular, studies on 
Tully-Fisher residuals applied to bright disk galaxies sug¬ 
gest that on average disk galaxies are sub-ma ximal with 
/c^(2.2 R^)=0.25-0.5 (Courteau & Rix 1999). Another 
approach assumes that the peak circular velocity of the 
stellar disk, measured at R=2.2 is related to the ver¬ 
tical velocity dispersio n, and the scale heig ht through 
some s caling relations ( van d er Kruit|1988). More re¬ 
cently, jBershady et al.| (|2010|) applied the velocity dis¬ 
persion tecEuIqueToTTsampe of 46 nearly face-on galax¬ 
ies finding that the disk fraction of their sample ranges 
between 0.16 and 0.5. A similar method has been re¬ 
cently applied to the Milky Way. The analysis includes 
the ability to disentangle stellar populations and results 
in a low estimate of the disk scale length and a corre¬ 
spond ingly high disk fraction of 0.69T0.07 (Bovy & Rix 


2013|) (hereafter BR13). 

The stud y here presents a model based on swing am¬ 
plification (Toomre |1981 ) to compute analytically the 
azimuthal wavenumber at each radius that is likely to 
experience the strongest growth to deter mine the num¬ 
ber o f spiral arms expected at that radiu s ( Athanassonla 
19881 [Tbomre fc Kalnais|[l991| |B osma|p^99| |4hchs &:| 


Athanassonla |2005p. For a sani ple of grand design spi- 


rals, Athanassonla et al. (1987), hereafter ABP87, ap 


plied the swing amplification formalism in an attempt 
to constrain the mass-to-light ratio of dozens of disks. 
ABP87 study is of considerable observational interest, 
since it treated the disk stars and gas as distinct compo¬ 
nents, and computed their contribution of the total ro¬ 
tation curve using directly their observed projected sur- 





















































2 


D’Onghia 


face density profiles with an appropriate potential solver. 
The bulge contribution was estimated in a similar way. 
This approach offered the advantage not to rely on a 
specific disk, bulge or halo model. It constrained the 
mass-to-light ratio of disk galaxies and thereby inferred 
the dark matter and disk contribution within the optical 
radius. Here the disk-stability criterion is generalized for 
a galaxy characterized by an exponential disk, a stellar 
bulge, and a Hernquist halo. It provides a simple ana¬ 
lytic expression to estimate the number of arms expected 
by swing amplification, making its application straight¬ 
forward to large surveys. 

The model is detailed in Section [2j Numerical simu¬ 
lations are presented in Section Section applies the 
method to an observational sample, and Section sum¬ 
marizes the main results. 


2. MODEL 

We consider the case of a stable disk with Q > 1 re¬ 
sponding to gravitational perturbations by swing amplifi¬ 
cation. For a given disk the efficiency of the amplification 
is characterized by a factor defined as: 
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Xcrit 27rGE m 
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where the critical wavelength is defined as 
^crit=47r^GE//^^, G is the gravitational constant, 
and K is the epicycle frequency. 

Numerical studies have shown that for a stable disk, 
as arms swing from leading to trailing, the amplification 
factor depends on the value of the X param eter. Partic¬ 
ularly suggestive is Fig. 8 of Toomre (1981) showing that 
for disks with flat circular velocity curves and Q = 1.2 
the gravitational response to perturbations of a stellar 
disk is the same as that of a gaseous disk and that the 
swing amplification is strongest when the parameter X 
is 1.5, l eading to an amplification of a f actor from 40 up 
to 100 (|Athanassoula||19841 |Fuchs||2001|) . 

Given the circular speed curve, the epicyclic frequency 
of the galaxy disk is determined. For a general galaxy 
with a bulge, disk and a dark halo, the epicyclic fre¬ 
quency depends on the total angular frequency, which is 
the sum of the angular frequencies, according to: 


= + (2) 

For an exponential disk, the angular frequency is: 
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where y = Rl2R(i^ Rd is the disk scale length, and Md 
is the disk total mass, including gas and stellar com¬ 
ponents. Fqr_^i^e^^;nd halos described by Hernquist 
models (Hernquist 1990), the angular frequency of the 
bulge is: 
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where and Mb are the bulge scale length and mass, 
respectively, and for the halo: 
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where and Mb are the halo scale length and the mass, 
respectively. 

The Fourier coefficient, m, that characterizes the num¬ 
ber of spiral arms, generalized for a galaxy with those 
properties becomes: 
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Fig. 1.— The total number of arms as a function of the distance 
from the center for a galaxy with the halo mass of 9.5x10^^ Mq 
and a halo scale length of 29.6 kpc. Galaxies with the same mass 
for the halo display a progressive increase of the disk fraction, 
from 20% to 50%, from the top to bottom. The predicted number 
of arms for t he Milky Way is also displayed adopting the structural 
properties of |Bovy fc Rix| ( |2013| > with the inclusion of a bulge with 
a mass of 4x10^ Mq and a scale radius of 600 pc. The assumed 
amplification factor is X=1.5. 


For a galaxy with a halo mass of 9.5x10^^ Mq the total 
number of arms m as a function of the distance from the 
galaxy center is displayed in Fig. 1. The different lines 
display models where the total mass distribution of the 
dark halo is kept fixed but the self-gravity of the disk is 
progressively increased. In particular, the disk mass frac¬ 
tion, fd^ within 2.2 scale lengths, rises from 20% to 50%, 
from the top to bottom. The disk scale lengths assumed 
are the same as in the numerical simulations presented 
in the next section. No bulge is assumed for simplicity 
except for the Milky Way model where parameters are 
assumed from BR13. 
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Note that the total number of arms is expected to in¬ 
crease with distance from the galaxy center as the highest 
Fourier components m, which are always present, domi¬ 
nate the spiral structures in the outer parts of the disk. 
This is in agreement with previous models and compar¬ 
ison with observations of ABP87. However, exponential 
disks are less self-gravitating as we move further out from 
the center, thus the arms in the outer parts of the disk 
will have lower strength a nd may just be in the form of 
poorly amplified wakelets (D’Onghia et al.||2013). 

Fig. 1 shows that at a given radius, the total number 
of spiral arms decreases as the disk fraction increases, as 
expected. The predicted number of arms of the Milky 
Way is also included in Fig. 1 (dashed line) adopting 
the structural parameters presented in BR13, with the 
disk fraction being estimated to be fd ?^0.69 within 2.2 
Rd. Note that according to this analysis the Milky Way 
should have two spiral arms at approximately 4.50 kpc, 
nearly at the edge of the stellar bar, and a total of 5-6 
spiral arms, lower in strength, in the solar neighborhood. 

3. NUMERICAL EXPERIMENTS 

Numerical simulations of various disk models have shown 
that changes to the disk mass lead to varying the crit¬ 
ical length scale parameter Xcrit lead to a differen t 
spiral morphology (see e.g. [S e llwood fc Car l bergI ( |1984|); 
Carlberg & Ereedman| (|1985P; INeison et ai.|(|2012Pb 
Here a set of numerical siniulations are enipioyed to de¬ 
termine the extent to which the response of the disk de¬ 
pends specifically on the mass distribution of the galaxy 
and on the disk self-gravity. 

In order to test the formalism introduced above, a set 
of numerical simulations are carried out with t he paral¬ 
l el Tre ePM code GADGET-3 (last descr ibed in Springel 
. Eollowing D’Onghia et al. (2013), the tree-based 


( 200 ^ ) 


gravity solver is employed coupled with a static poten¬ 
tial to solve for the evolution of collisionless particles. 
The galaxies in this study consist of dark matter ha¬ 
los and rotationally supported stellar disks of 10^ par¬ 
ticles. The parameters describing each component are 
independent and the mode ls are const r ucted similarly to 
the approach described in Hernquist (1993). The dark 
matter mass distribution is modeled assuming a Hern¬ 
quist profile with a total halo mass and a scale length of 
9.5x10^^ Mq and 29.6 kpc, respectively. The stellar disk 
is modeled according to the initial conditions as a thin 
exponential surface density profile of scale length so 
that the stellar disk mass is = fd^tot where fd is the 
disk fraction of the total mass of the galaxy. The vertical 
mass distribution of the stars in the disk is specified by 
the profile of an isothermal sheet with a radially constant 
scale height zq. A value of = 0.1 R^^ is adopted. 

In a sequence of experiments shown in Eig. 2, the 
contribution of the disk to the total mass of the galaxy 
is altered. Like in Eig. 1 the total mass of the halo is 
fixed but the disk mass is increased in order to obtain 
a disk fraction, /^, within 2.2 disk scale lengths varying 
from 20% to 50% of the total mass. Disk scale lengths 
are set in the initial conditions to values that conserve 
the total angular momentum of the disk. They range 
from 2.75 kpc in the model with fd=0.5 to 3.3 kpc in 
the one with fd=0.2. Eig.2 shows that the galaxy ranges 
from the multi-armed spirals obtained for disks with low 
mass fractions (lower panels in Eig. 2), to those with 



Fig. 2. — Spiral morphologies of the simulated stellar disk dis¬ 
played face-on after one galactic year. The disk mass fraction at 
2.2 scale lengths, fd ranges from 50% of the total mass (top left 
panel) to 20% (bottom right panel). 


four prominent arms typical of galaxies where the disk 
contributes about 50% of the total mass (top left panel 
in Eig. 2). As anticipated by the analytic arguments and 
previous studies, the results shown in Eig. 2 confirm that 
the spiral morphology in the simulated disks is indeed 
determined by the structural properties of the galaxy 
through the expression of Xcrit- 

A quantitative analysis of the strength of the arms as a 
function of the galaxy radius has been performed through 
a Eourier analysis of the surface mass density of the stel¬ 
lar disk. The results show that the amplitude of the 
arms, on average, ranges from 10% of the stellar back¬ 
ground for the multi-armed galaxies (right bottom panel 
of Eig. 2) to 50% for the galaxies with four arms and ap¬ 
proximately four-fold rotational symmetry (top left panel 
in Eig. 2). The number of arms increases as a function 
of radius in agreement with the analytic estimates. 


4. APPLICATION TO THE DISKMASS SAMPLE 

The formalism introduced can then be applied to the 
Milky Way and to the DiskMass sample, consisting of 
face-on galaxies with rotation velocities between 100 
kms“^ and 250 kms“^ ( Bershady et al.|2010 ; Davis et al. 
2015). Models are run assuming the structural proper¬ 


ties of the galaxies, when available, i n particular scale 
length s, and total mass as reported by [Martinsson et ah 
(2013), with the observed disk fraction fd ranging from 
0.16 to 0.5. As expected, Eig. 3 shows a correlation be¬ 
tween the total number of spiral arms predicted by the 
swing amplification and the disk fraction. Submaximal 
disks are thus expected to be multi-armed galaxies, and 
the two-armed galaxies are the ones with higher fd within 
2.2 scale lengths. A visual inspection of the images of the 
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Fig. 3. — Total number of spiral arms predicted by models ap¬ 
plied to the DiskMass catalog of galaxies (blue symbols) as a func¬ 
tion of the disk fraction, measured at 2.2 disk scale lengths, m is 
the average value estimated by models assuming 1.5 < X <2 with 
the error bars given by the two extreme values. The prediction of 
the number of spiral arms of the Milky Way is also included (green 
sy mbol) for the struc tural parameters of the Milky Way reported 
in |Bovy &: Rix| ( |2013| >. 


galaxy sample in the available bands supports this result, 
although a more quantitative analysis is needed. 

5. CONCLUSIONS 

A model based on the swing amplification mechanism 
is generalized for a realistic galaxy model with an ex¬ 
ponential stellar disk and a dark halo described by the 
Hernquist mass profile. The number of spiral arms is 
derived in the context of disk stability. More generally, 
the model confirms previous results that show that the 
number of arms in galaxies should depend only on the 
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structural properties of the galaxy, with a correlation 
between the disk fraction and the dominant wave mode 
numbers. These models are aimed at providing an ad¬ 
ditional constraint on the number of spiral arms in the 
Milky Way as a function of the radius and an indepen¬ 
dent test on the structural parameters of the Milky Way 
provided by BR13 (as shown in Fig.3). They have been 
applied to the DiskMass survey to confirm the known 
correlation between the total number of spiral arms pre¬ 
dicted by the swing amplification theory and the disk 
mass fraction. 

Ongoing and forthcoming surveys of disk galaxies pro¬ 
vide strong tests of the applicability of these models. In 
particular, they can constrain the mass-to-light ratio of 
disk galaxies and hence to infer the dark matter and 
disk contribution within the optical radius. Although the 
study presented in ABP87 is model independent and with 
a more accurate treatment of the swing amplification, the 
formalism presented here has the advantage to provide an 
analytic expression for the number of arms expected by 
swing amplification that can be used to compare with the 
actual number of arms observed in large sample of galax¬ 
ies. For a given rotation curve the epicycle frequency is 
determined analytically. If the number of arms expressed 
by m at a given radius can be measured independently, 
that number indicates the azimuthal wavenumber un¬ 
dergoing the strongest swing amplification at that radius 
(according to eq.(6) for X = 1.5 —2 for any given galaxy). 
Therefore the surface density of the disk can be derived 
analytically from eq.(l), leading to an independent con¬ 
straint in the mass-to-light ratio of the disk, and hence 
to the halo contribution to the total mass. 
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